Vortex Ring Reconnections 
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We investigate numerically the Navier-Stokes dynamics of reconnecting vortex rings at small Re 
number. We find that reconnections are dissipative due to the smoothing of vorticity gradients at 
reconnection kinks and to the formation of secondary structures of stretched anti-parallel vorticity 
which transfer kinetic energy to small scales where it is subsequently dissipated efficiently. In 
addition, the relaxation of the reconnection kinks excites Kelvin waves which due to strong damping 
are of low wavenumber and affect directly only large scale properties of the flow. 

PACS numbers: 47.32.Cc, 02.70.Ns 
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In flow phenomena as diverse as quantum Q magnetic 
H and incompressible [[§ fluids, it is useful to study the 
physics of turbulence by modeling the system as a col- 
lection of tubular flux loops which in the case of vortical 
fields are called vortex filaments. An intrinsic property 
of such highly structured systems is their ability to dy- 
namically change their topology via reconnection mecha- 
nisms. Does this change in topology affect in turn prop- 
erties of fluid turbulence like the intermittcncy and the 
scalar-mixing (which depend directly on the structure of 
the flow) or like the dynamics of energy in wavenum- 
ber space? Or is it the case that reconnection events 
are not generic and thus have no direct impact on the 
mean properties of turbulent flows? The aim of this letter 
is to address these issues by fully resolving the Navier- 
Stokes dynamics of interacting vortex rings for three sim- 
ple geometries having great potential for illuminating the 
physics of reconnection. Although the flows considered 
are not strictly turbulent, the hope is that in a future 
structural approach to the problem of turbulence a sig- 
nificant part of the flow complexity could be traced back 
to the physics of similar vortex interactions. 
Incompressible vortex reconnections have an extensive 
bibliography (for a review of the work up to 1994, see 
In (|, || reconnections of vortex tubes were consid- 
ered with an emphasis on the possibility of singularity 
formation as Re — * oo. In JtJ the strong interactions 
between vortex rings were computed with the interest 
in developing numerical methods and turbulence models 
rather than in focusing on the physics of reconnection. 
In H it is discussed how a linked vortex configuration 
could be achieved starting from an unlinked initial state 
and in it is considered how the mixing of a nondiffus- 
ing passive scalar is affected during vortex ring collision. 
The reconnection of two approaching (but not colliding) 
vortex rings was studied experimentally in ]10| | and the- 
oretically in p] ]. This letter extends these studies by 
considering generic vortex configurations and by captur- 
ing more features of vortex reconnections in a turbulent 
flow. 



We solve the Navier-Stokes equations for an un- 
bounded three-dimensional incompressible viscous flow. 
We employ the vorticity formulation: 



u ■ V ) w = (Vu) • uj + vV 2 ui , (1) 
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V-it = 0, (2) 

where u is the velocity and u> is the vorticity. We use a 
vortex particle method |L2j . In this method, the vorticity 
is discretized with Lagrangian elements. These elements 
which carry a vector- valued Gaussian distribution of vor- 
ticity are convected and stretched by the local velocity 
obtained by the Biot-Savart law. The complexity of the 
velocity computation is normally 0(N 2 ) with N being 
the number of particles; we have used a multipole algo- 
rithm that reduces this complexity to 0(N\og(N)). The 
viscous diffusion is modeled by the Particle Strength Ex- 
change scheme. 

We calculate some global quantities: the kinetic energy 
E, the enstrophy fl and the helicity H , defined as 



E = - I u ■ udx . 
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£1 = / u> ■ u> dx , 



(3) 
(4) 
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For unbounded flows, the relation between kinetic energy 
and enstrophy is 



H = I uj ■ udx . 



dt 
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(6) 



We also compute the evolution of the spectrum of the 
kinetic energy E(k) which, in terms of the Fourier trans- 
form of vorticity uj — J u>(r) e ~ lr k dr, is defined 



E{k) = -{2kY I uj-u*dn k 



(7) 
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where dO^ denotes sin 9^ d9k d<pk , the solid angle element 
in spherical coordinates. The calculation of the spectrum 
requires a double summation over the vortex elements 
which results to 0(N 2 ) complexity. Because of this, the 
calculation of the spectrum is much more costly than 
the solution of the Biot-Savart law. Since the number 
of particles grows substantially during our simulations, 
from around N = 5 10 4 at t = to 8 10 5 in the end, our 
computational resources did not allow us to compute the 
spectra for all times. 

We consider three reconnection configurations at the 
same Reynolds number: Re = E = 250 where T is the 
circulation of one ring and v is the kinematic viscosity. 
This small value of the Re was dictated by the compu- 
tational cost and the need for well-resolved reconnection 
regions. All the rings have the same initial T. All of our 
conclusions are conditioned upon the small value of the 
Re number, as well as, on the common initial circulation 
and should not be extrapolated uncritically to other set- 
tings. The vorticity distribution in the cross-section of 
every ring is Gaussian with a cut-off 

where r is the distance to the core center, a is the 
core radius and ojg is the azimuthal vorticity. We chose 
a = 0.05 R (where R is the radius of the ring), to en- 
sure that the rings are still thin when reconnections oc- 
cur. Our results were made dimensionless in the following 
manner: t — x = ^, uj = R ^ where t', x', w' are 
dimensional. 

We study three configurations. In the first case (Fig. 
[I], H and ||) , the initial rings are placed at a distance of 
i?/4 apart in the z direction, offset by R along the y axis 
and they move in opposite directions along the z axis. 
The second geometry (Fig. |[ || and ||), consists of two 
rings of different radii (i? and R/2) moving in the same 
direction along the z axis, with the center of the small 




FIG. 3: Vortex rings in an offset collision: evolution of the 
energy spectrum 

ring in collision course with the circumference of the large 
one. The small ring has a larger self-induced velocity and 
catches up with the large ring. Finally, in case three (Fig. 
@, | and ||), the two rings are linked at 90° a ring going 
through the other in its center. One is moving in the pos- 
itive z direction; the other, in the positive y direction. 

All three evolutions have basic common features which 
will be discussed in the context of the first configura- 
tion for clarity. The spectrum at t = (Fig.||) has the 
characteristic oscillations of the spectrum of isolated vor- 
tex rings and a cuf-off at the scale of ring core radius 
a = 0.05, k = 20. Our data suggest that the reconnec- 
tion starts around t = 0.6 when an acceleration of energy 
decay appears, and ends around t = 1.4. Specifically, as 
the rings approach each other, they stretch and deform 
near the collision points so that their respective vortic- 
ities become locally anti-parallel. The two ends of this 
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FIG. 4: Vortex rings of different radii: contours of vorticity; 

, the contour is lo = 0.15 6 -t ~ 
it is w = 0.05 w*=° 



FIG. 5: Vortex rings of different radii: kinetic energy and 
enstrophy 




FIG. 6: Vortex rings of different radii: evolution of the energy 
spectrum 



stretching region eventually become reconnection kinks 
in which (in the absence of singularities), the strong vor- 
ticity gradients are smoothed out by diffusion. This is 
also seen in the graphs of the global quantities (FigJ|) 
where the beginning of the reconnection process corre- 
sponds to a hump in the graph of and to a steepening 
of the slope of E(t) between t = 0.6 and t = 1.4. The time 
length of the reconnection is significantly longer than the 
viscous time scale t v — a 2 jv — 6.25 10~ 2 and of the same 
order as the convective time t c = R 2 /Y — 1. This con- 
tradicts the findings of |l(| where it is reported that the 
viscous scale is much larger than the reconnection dura- 
tion. However, in [|l0| the Re number was 1600. The 
conclusion that the duration of reconnection is inversely 
proportional to the Re number and thus to the circula- 
tion of the vortices is plausible (also in agreement with 
H ^3|), but it is subject to the condition in that the 



rings are merely touching themselves rather than collid- 
ing. 

After some time (t = 1.6), we can say that two new rings 
are formed. The pairs of filaments between the reconnec- 
tion regions are stretched further as the new rings move 
apart from each other (t = 1.6 to 5.6). Viscous diffusion 
weakens the vorticity magnitude in these structures while 
the reconnection kinks relax in the form of low wavenum- 
ber Kelvin waves. These stretched vorticity structures 
are responsible for a continued transfer of energy to the 
smallest scales until these structures are dissipated away. 
This we conclude by noticing that the high wavenumber 
cut-off of the spectrum becomes a non-exponential one 
(although still very steep) and that between t = 0.8 and 
2.4 there is a significant decay of the energy spectrum for 
k < 20 but little change for k > 20. 

It is conceivable that for Re numbers higher than 250 an 
intermediate scaling range (that is in between the k~ l 
and k~ 7 regimes) could appear with inertial type of scal- 
ing. It is also expected that with increasing Re number 
the hump observed in the global enstrophy during the re- 
connection could become more pronounced and (accord- 
ing to the previous discussion) shorter in duration. 
The initial helicity is non-zero in the linked configura- 
tion. As the rings evolve and unlink themselves, the cor- 
responding helicity increases as expected and one clearly 
notices a slope steepening between t—3 and 5 (Fig.||) 
when the reconnection action occurs in this case. 

In conclusion, we studied three generic vortex ring 
configurations and we found that in all cases the rings 
reconnect. This suggests that reconnection is a common 
phenomenon in vortex filament encounters and perhaps 
also in turbulent flows. In addition, we observe an 
intensification of dissipation which is local in time 
and could be a mechanism contributing to turbulence 
intermittency. A by-product of reconnection is the 
formation of stretched structures with anti-parallel 
vorticity which transfer energy to the smallest scales 
where it is rapidly dissipated. Without this energy 
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FIG. 9: Linked vortex rings: evolution of the energy spectrum 

redistribution in wavenumber space the decay of global 
kinetic energy would have been slower. This important 
effect depends directly on the details of the initial vortex 
configuration (compare with experiments in |lC[|). The 
observed intensification of small scale motions hints to 
an enhancement of small scale mixing of passive scalars 
with Sc > 1. The excited Kelvin waves represent a fast 
mechanism for energy transfer, but the small Re number 
of our calculations is not suitable for understanding 
their full importance. In particular, they are confined 
to low wavenumbers in opposition to the Kelvin waves 
observed in reconnections in quantum fluids JlJ] . This is 
because quantum filaments are inviscid and have a very 
thin core (a ~ 0.1 nm) so that high wavenumber Kelvin 
waves propagate without damping even for rings with 
small circulation. 



Besides illuminating important physics, the present 
work will guide future introduction of phenomenological 
reconnection models into vortex filament computational 
methods. In this way, the applicability of the latter 
methods will be extended to flows with complex vorticity 
configurations. 
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